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Fault location using measurements of current and voltage from 
5 one end of a line 

TECHNICAL FIELD. 

The present invention is concerned with a fault location 
technique for a section of a power transmission line utilizing 
10 measurements of current and voltage made at terminals located 
at one end of the section of the power line. 



BACKGROUND ART 

Several methods and approaches for fault location in high 
15 voltage power transmission systems, and power distribution 
systems, have been developed and employed. One approach has 
been to use voltage/current transducers located at terminals 
located at each of two ends of a section of the power line to 
be monitored. Inductive current transformers are used to 
20 provide a measurement of instantaneous current in a 
transmission line. 

US 4,559,491 entitled Method and device for locating a fault 
point on a three-phase transmission line, describes a one-end 

25 fault location (FL) algorithm. High accuracy of fault location 
using a fault locator device at one end of a line is achieved 
by taking into account the actual distribution of a fault 
current in the transmission network. This algorithm has been 
successfully implemented into a product in 1982 and is in 

30 operation with single and parallel transmission lines in many 
countries around the world. However, for certain fault 
conditions it is difficult to obtain accurate pre-fault 
quantities, such as pre-fault currents, in order to calculate 
an estimate for voltage drop across the fault path. Also, a 



disadvantage of using phase voltages and currents and zero 
sequence components of currents is that it is relatively 
difficult using these vaiues to compensate for shunt 
capacitance effects. In addition, the fault locator method 
5 described is not suitable for single and parallel line 
sections which have an extra link across the ends of the 
sections . 

SUMMARY OF THE INVENTION 
10 The aim of the present invention is to remedy one or more of 
the above mentioned problems. 

This is obtained by a method characterised by claim 1. 
Specific features of the present invention are characterised 
15 by the appending claims. 

In one aspect of the invention, a method comprising a new 
formulation of a one-end fault locator algorithm has been 
proposed. The uniform description of the transmission network 
20 in terms of symmetrical components as well as the generalized 
models of fault loops and faults have been applied. The 
resulting advantages include, the algorithm can be uped for 
locating faults in typical single and parallel transmission 
lines, and, in addition, fault location may also be carried 
25 out for both single and parallel lines with an extra link 

between the line ends. Another advantage is that a procedure 
: for calculation of a distance to fault is in the form of a 

compact quadratic equation with the coefficients dependent on 
: a fault type, acquired measurements and impedance data for the 

:30 transmission network. Another advantage of the invention is 
: that optimal estimation of the voltage drop across a fault 

: path is applied, which has the result that the pre-fault 

: currents in case of single phase-to-ground faults and phase- 

to-phase faults are no longer required. 
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In an embodiment, compensation for shunt capacitances is 
facilitated by means of the use of the notation of symmetrical 
components. The distributed long line model of the line has 
been applied for that, The compensation is performed 
5 individually for all the sequences . The currents for 

particular sequences are compensated for the shunt currents 
and then the fault loop compensated current is composed. In 
another embodiment improved accuracy has been obtained by 
means of sin option to measure the source impedance at the 
10 remote end instead of using a representative value. The source 
impedance measured at the remote end may be considered as sent 
to the fault locator by using a simple communication means. 

In another embodiment, a method for one end fault location for 
parallel lines to locate single phase- to-ground faults is 
15 described under a plurality of conditions. In another further 
embodiment a method is described for one end fault location 
with standard availability of the measured signals for ground 
faults including both single phase-to-ground faults and phase- 
to-phase- to-ground faults. 

20 

In , another, aspect of the invention, a fault locator device for 
carrying out the method of the invention is characterised by 
claim 18. Specific features of the inventive fault locator 
device are characterised by the appending claims . 

25 

In another aspect of the invention a computer program is 
described for carrying out the method according to the 
invention. In another aspect of the invention a computer 
program product comprising a computer program for carrying out 
3 0 the method of the invention is described. In another, further 
aspect of the invention a graphical user interface is 
described for displaying a distance to a fault from one end of 
a section of a power line. 



BRIEF DESCRIPTION OF THE DRAWINGS 



A more complete understanding of the method and system of the 
present invention may be had by reference to the following 
5 detailed description when taken in conjunction with the 
accompanying drawings wherein: 

Figure 1 shows in a single schematic diagram a method of fault 
location in power transmission and/or distribution systems for 
10 parallel lines and single lines according to an embodiment of 
the invention; 

Figure 2a shows a schematic circuit diagram for a parallel 
transmission network for the positive sequence component in 
15 which the fault loop is marked for the case of a fault locator 
installed at the terminal AA. Figures 2b, 2c show 
corresponding diagrams for the negative sequence and zero 
sequence components, respectively; 

20 Figure 3a is a schematic block diagram for obtaining and 
calculating the phasors of the symmetrical components of 
voltages and currents used for composing the fault loop 
voltage. Figure 3b shows a corresponding diagram for composing 
the fault loop current; 

25 

Figure 4 shows a circuit diagram for determining the fault 
current distribution factor for the positive sequence of a 
single line, in which diagram quantities for the negative 
sequence are shown indicated in brackets; 

30 

Figure 5 shows a circuit diagram corresponding to Figure 4 for 
single lines for determining the fault current distribution 
factor for the positive sequence of parallel lines, in which 
quantities for the negative sequence are also shown indicated 
35 in brackets; 
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Figure 6 shows a schematic diagram for an embodiment of the 
invention in which source impedance measured at a remote end B 
may be communicated to a fault locator at the first end A; 



5 Figure 7 is a circuit diagram of an embodiment in which the 
shunt capacitances are taken into account, and shows a 
positive sequence circuit diagram during a first iteration; 

Figure 8 shows a negative sequence circuit diagram for taking 
10 the shunt capacitances effect into account during a first 
iteration; 

Figure 9 shows a zero sequence circuit diagram for taking the 
shunt capacitances effect into account during a first 
15 iteration; 

Figure 10 shows a flowchart for a method for locating a fault 
in a single line according to an embodiment of the invention; 

20 Figure 11 shows a flowchart for a method for locating a fault 
in parallel lines according to an embodiment of the invention; 

Figure 12 and Figures 13a, 13b, 14, 15a and 15b show schematic 
diagrams of possible fault- types (phase-to-phase, phase- to- 
25 ground and so on) with respect to derivation of coefficients 

* r r 

: *; : for Table 2 in Appendix A2 . Figure 12 shows fault types from 

: 4 a-g, and Figures 13a, 13b faults between phases a-b. Figure 14 

[ shows an a-b-g fault. Figures 15a and 15b show symmetrical 

faults a-b-c and a-b-c-g respectively; 

* r m 

* * * 

• » 

... Figures 16 and 17 show schematic diagrams for the derivation 

l[\ of the complex coefficients in the fault current distribution 

» » * 

r./ factors for the positive (negative) sequence included in Table 

* ■ 

3. Figure 16 shows the case of a single line with an extra 



link between the substations . Figure 17 shows the case of 
parallel lines with an extra link between the substations; 

Figure 18 shows a fault locator device and system according to 
an embodiment of the invention; 

Figure 19 shows a flowchart for a method for locating a single 
phase -to -ground fault in parallel lines in the case of 
measurements from the healthy line being unavailable, 
according to an embodiment of the invention; 

Figure 20 shows a schematic diagram for a method of fault 
location for parallel lines with different modes of the 
healthy parallel operation; 

Figures 21a, shows a schematic equivalent circuit diagram for 
a parallel network for the incremental positive or the 
negative sequence. Figure 21b shows an equivalent circuit 
diagram for the zero sequence while both parallel lines are in 
operation. Figure 21c shows the equivalent circuit diagram for 
the zero sequence with the healthy parallel line switched off 
and grounded; 

Figure 22 shows a flowchart for a method for locating phase- 
to-phase and phase-to-ground faults in parallel lines in the 
case of providing the zero sequence currents from the healthy 
parallel line according to another embodiment of the 
invention; 

Figure 23 shows a schematic diagram fault location for 
parallel lines with standard availability of measurements 
according to another embodiment of the invention; 



Figures 24 a,b,c show the equivalent circuit diagrams of 
parallel lines for positive, negative and zero sequence 
currents respectively . 



DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Fig.l presents a schematic diagram for one-end fault location 
applied for parallel lines and for a power transmission or 
distribution system with a single line. A fault locator 1 is 
positioned at one end 2 of a single line AA-BA 3 or parallel 
lines AA-BA, AB-BB, 4. A fault F is shown at FA with a 
corresponding fault resistance, 5, denoted as Rf . A value for 
distance to the fault d from one end 2 determined and provided 
by the fault locator 1 is indicated with the reference number 
7. Components such as parallel line AB-BB and quantities such 
as a parallel line value zero sequence current Iabo shown with 
dotted lines are excluded when considering a single line case. 

The fault locator 1 positioned at the first end 2, or *A' end, 
is supplied with the following input signals: 

- three-phase voltages (Vaa) of the faulted line, 

- three-phase voltages (Iaa) of the faulted line, 

- zero sequence current (Iabo) % from the healthy parallel line 
(zero sequence is not present when the single line only is 
considered) . 

Fig. 2 a,b,c show circuit diagrams of a parallel transmission 
network for the positive 2a, negative 2b, and zero sequence 2c 
components. Fault loops for the sequence components 21a, 21b, 
21c are shown for the case of a fault locator installed at the 
terminal AA. An extra link 25 between the terminals A, B is 
shown. A generalized model of the fault loop considered for 
different fault types is stated as: 

V Mj -dZ l£A LMj> - r f(&fiL F] + « F1 l F2 +£foIfo) = 0 < 1} 
where : 
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d - unknown and sought distance to fault, 
Z lIA ~ positive sequence impedance of the faulted line, 
V I - fault loop voltage and current composed according 

— A4 _p ' —AA_p 

to the fault type, 

R F - fault resistance, 

l F . - sequence components of the total fault current (i = 0, i 
=1, i = 2) , 

a Fi - weighting coefficients (TABLE 2) . 

Fault loop voltage and current can be expressed in terms of 
the symmetrical components of measured voltages/ currents : 
V MJS> =a,V Ml +a 2 V M2 +a 0 V MO {2) 

z z 



(3) 



where : 

L5 AA, AB - subscripts used for indicating measurements 

acquired from the faulted line (AA) and from the healthy line 
(AB) , respectively, 

a*. *x> 2-1 ~ coefficients which are gathered in TABLE 1 (the 
Tables are arranged below at the end of the description of 
embodiments and derivation of these coefficients is shown in 
APPENDIX Al, also attached) . 

Zola > Zom - impedance of the faulted line and mutual coupling 
between the lines for the zero sequence, respectively, 
a 0m - o.o ~ for parallel lines, 
25 a 0m =0 - for single lines. 



The phasors of symmetrical components of voltages, positive: 
V^,, negative: V MZ and zero sequence: V M0 as well as the 
phasors of symmetrical components of currents, positive: Laa\> 
30 negative: / M2 , zero sequence from the faulted line: Lmo a nd 
zero sequence from the healthy line: l AB0 are calculated from 



the acquired measurements as shown schematically in schematic 
block diagrams Figures 3a and 3b. 



Figure 3a shows an input of instantaneous phase voltages 3 0a, 
5 filtering stage 33a, phasors of phase voltages 31a, 

calculation of phasors of symmetrical components 33b and 
phasors of symmetrical components of voltages output at 32a. 
It may be seen from Fig 3a that acquired phase voltage 
measurements are subjected to a filter, then calculations are 

10 made to find the symmetrical components of the fault loop 

voltage. Figure 3b shows correspondingly stages used to find 
the symmetrical components of the fault loop current. Figure 
3b shows instantaneous phase currents and instantaneous zero 
sequence current from the healthy line 30b, filtering 33b, 

15 phasors of phase currents and phasor of zero sequence current 
from the healthy line 31b, calculation 34b and phasors of 
symmetrical components of currents output at 32b. 

Fault loop signals may be composed according to formulae (2)- 
20 (3) and TABLE 1, which is the alternative to the classic 

approach (TABLE 1A, fault loop voltage (Zaa_ p > and current 
(l M p ), which was used in the fault locator from [1-2]. 

Voltage drop across a fault path resistance, the third term in 
25 (1), can be expressed in terms of the current distribution 

factors and local measurements of currents which results in: 

Formula (4) has been obtained from the following relations 
3 0 between the symmetrical components of a total fault current 
and measured currents: 






(5) 



where : 



l Fi ; / F2 ; l F0 - symmetrical components of a total fault 
5 current , 

k Fl , k F2 , k FQ - fault current distribution factors for 
particular sequence quantities, 

AlAAi=I/Wi-We; LaaH Laao ' symmetrical components of 
currents measured in the line A at the station A (subscript 
10 AA) ; note that in case of the positive sequence the 

incremental quantity (post-fault current minus pre-fault 
current) is used. 

Voltage drop across the fault path, as shown in the third term 
15 in equation (1), is expressed using sequence components of the 
total fault current. The weighting coefficients a F0 , a Fl , a F2 , 
can accordingly be determined by taking the boundary 
conditions for particular fault type. See TABLE 2, Alternative 
sets of the weighting coefficients for determining a voltage 
20 drop across the fault path resistance. Examples of derivation 
of these coefficients are contained in APPENDIX A2 . 

There is some freedom for setting the weighting coefficients. 

It is proposed to utilize this freedom firstly for avoiding 
: 25 zero sequence quantities, since the zero sequence impedance of 

a line may be considered as an unreliable parameter. This can 
\ be accomplished by setting a F0 =0 as shown in TABLE 2. 

: Secondly, the freedom in establishing the weighting 

. 3 0 coefficients can be utilized for determining the preference 

: for using particular quantities. The negative sequence { TABLE 

: 2, set I) or the positive sequence (TABLE 2, set II) can be 




preferred as well as possibly both types of the quantities 
(TABLE 2, set III) can be used for determining the voltage 
drop across the fault path. 

5 The set I is recommended for further use, thus avoiding the 
positive sequence, and thus avoiding the pre-fault positive 
sequence current, for the largest number of faults. Avoiding 
the pre-fault positive sequence current is highly desirable 
since sometimes the pre-fault currents - due to certain 

10 reasons - can not be recorded or registered, but may be 

contaminated by one or more the symptoms of the occurring 
fault. Moreover, the accuracy of recording the pre-fault 
currents, which are basically lower than the post-fault 
currents, is not very great. This is so since the A/D 

15 converters operate with less accuracy in the low range. 

Fault current distribution factors depend on the configuration 
of the transmission network, Figures 4, 5, and impedance 
parameters. Basically, all impedances for the positive and for 
20 the negative sequence are equal to each other and thus one 
obtains : 

k £ ■= (6) 

Coefficients in a fault current distribution factor (6) for a 
single (Figure 4) and for parallel lines (Figure 5) are 
25 gathered in TABLE 3, Coefficients for determining a fault 
current distribution factor, (note that derivation of the 
coefficients is shown in APPENDIX A3). 



Figure 4 shows a circuit diagram of a single line for 
30 determining the fault current distribution factor 

for the positive sequence currents and with the negative 
sequence currents such as shown in brackets . Similarly 
Figure 5 shows a circuit diagram of parallel lines for 
determining the fault current distribution factor with 
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positive sequence currents wherein the negative sequence 
currents are shown in brackets. 

In Fig. 4 the extra link 45 between the terminals A, B having 
5 impedance for the positive sequence equal to: Z lAB can be 

considered as existing {Z lAB *~ ) or as not present (Z lAB ^>~) . 
in Fig. 5 the extra link 55 between the terminals A, B having 
impedance for the positive sequence equal to: Z lAB can be 

considered as existing ( Z lLB&AB = fr^ 5 ) or as not present 

Substituting (6) into (4) and adjusting £ F0 =0 (as in TABLE 2) 
results in: 

After multiplying both sides of (7) by: K f +Li and some 

LAA_p 

15 rearrangements, the quadratic formula with two unknowns, d- 

[p.u.] sought fault distance from A, R F - fault resistance, is 
obtained : 

—AAjp 



20 



where : 



Zaa = — dd=B - calculated fault loop impedance. 



Writing formula (8) in more compact form results in: 

A 2 d 2 +A l d + A 0 + A QO R F =0 {8a) 

where : 

25 A 2 = A 2Re ^]A 2lm =K x Z llA 

A\ = A l Kc +jA t _ lm ~L\Z\LA -K\Z.AA_p 
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Z,, ^rr^^S - calculated fault loop impedance 
£i * i M.i " coefficients gathered in TABLE 3 . 

5 Formula (8a) can be written separately for real and imaginary 
parts : 

^2_Re<* 2 + A_Re<* + A>_Re + 4)0_R«*F =0 < 8b) 

A^d 2 + A llm d + V.lm + A)0„lm^ =0 ( 8C > 

10 Combining (8b) and (8c) in such the way that fault resistance 
is eliminated [that is, equation (8b) is multiplied by i4o 0 _, m 
and equation (8c) by 4)o_Re and then subtracting them] yields 
the quadratic formula for a sought fault distance: 
B 2 d 2 +B l d+B 0 =0 (9) 

15 where : 

B\ — A_ReA)0_Im ~ A„ImA)0_Re 
#0 = A)__ReA)0_Im ~ A)_JmA)0_Re 

20 Equation (9) has two roots (dj, d 2 ) for the distance to fault: 



— J5 t — ^Bf —4B 2 Bq 



2# 



2 



. -B] + ^B?-4 B 2 B 0 

do = 



2 2S 2 



(10) 



The root which fulfils the condition (0<d<l) is selected as 
25 the solution for the distance to fault. 
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In another embodiment of the invention, the method of fault 
location is carried out by using a measurement of source 
impedance at the second end remote from the fault locator 1. 
instead of a representative value for the source impedance at 

5 the remote end, and by communicating that measured value to 
the local end using a communication means. Coefficients from 
(9) are determined with the local measurements and the 
impedance data for the transmission line, the extra link 
between the line terminals and the equivalent systems at the 
10 line terminals. Impedance of the equivalent system at the 
local substation (Z XsA ) can be traced on-line with the local 
measurements. In contrast, the remote system impedance (Z ljB > 
is not measurable locally from A. Thus, the "representative" 
value of this impedance may be provided for the algorithm [1- 

15 2] . 

The alternative solution for the single line case is shown in 
Figure 6, which shows a fault locator 1 at the first end 2 
near to system A, and another device 10 located at the remote 
20 end 3 close to system B, indicated as RD. A communication 

signal 9 is shown being sent from the 10 at the remote end to 
the fault locator 1 at the local end. 

The remote source impedance (Z ]sB ) is measured by the remote 
25 device RD, 10, which may be another fault locator or any 
suitable device such as a digital relay or digital fault 
: recorder, in the remote substation and the measurement 9 is 
sent via a communication channel 60. Synchronization of 
measurements at the line terminals is not required. The source 
3 0 impedance is calculated from the known relation between the 
incremental positive voltage lAV m ) and the incremental 
positive sequence current (4/ B i) [3-4]: 

7 _ AZbi (ID 
- UB ~ AX 



15 

Similarly, the fault locator 1 calculates the local source 
impedance 



5 In another and preferred embodiment of the invention , 

compensation is carried out for the shunt capacitance of the 
line. Compensation for shunt capacitances effect can be 
accomplished by taking into account the lumped model (only the 
longitudinal R-X L parameters are taken into account) or the 
10 distributed long transmission line model- The distributed long 
line model [5] as providing higher accuracy of fault location, 
has been considered here. 



The compensation for the single line is presented further. 
15 This means that when composing fault loop current (3) the term 
reflecting the mutual coupling effect disappears (a 0/n -0). 
Moreover, the single subscript (A instead of AA) is used. 



Fault location procedure with compensating for shunt 
20 capacitances of a transmission line requires the following 
additional input data: 

Cj£ - shunt capacitance of a whole line for the positive and 
the negative sequences (parameters of a line for the positive 
and the negative sequences are identical and thus; C 2 i,=C JL ) 
25 c ol " shunt capacitance of a whole line for the zero sequence, 
I - total line length (km) , used for expressing 
impedances /capacitances of the line per km length. 

The compensation of shunt capacitances may be introduced while 
3 0 determining the voltage drop across the faulted line segment - 
the second term in the generalized fault loop model (1) . This 
requires compensating the components of the calculated 
currents for particular sequences. Thus, the original measured 
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currents : Lai > Lai > L A o have to be replaced by the currents 
after the introduced compensation: LA\_comp > Lhi_camp • LAO_a»np • 
At the same time the original fault loop voltage, the first 
term in the model (1), is taken for a distance to fault 
5 calculation. As concerns determining the voltage drop across 
the fault resistance, the third term in (1) , it is assumed 
here, which is a standard practice, that the effect of line 
capacitances at the location of the fault (point F) , may be 
neglected. This is justified as the impedance of the 
10 capacitive branch at that location is much greater than the 

fault resistance. This means that the voltage drop across the 
fault resistance is determined without taking into account the 
shunt capacitances . . . 

15 Calculating a distance to fault the following impedances 
(defined below) are taken as: 

Z{T* - positive sequence impedance of a line with taking into 
account the distributed long line model, 
Z$L ~ as above « but for the zero sequence. 

20 

The.. compensation procedure requires iterative calculations, 
performed until the convergence is achieved (i.e. iterated 
until the location estimate ceases to differ from the previous 
estimate) . However, the studies conducted revealed that 
25 results of acceptable accuracy may be obtained using a fixed 
: number of iterations, for example, 2-3 iterations. The 

calculated distance to a fault from a particular (say, present 
iteration) is utilized for determining the shunt current in 
the next iteration. The determined shunt current is then 
: 30 deduced from the measured current. A distance to fault, 

calculated without considering the shunt effect (10) , is taken 
as the starting value for the first iteration. 




A way of conducting the first iteration of the compensation is 
shown in Figures 7, 8, 9 for the positive sequence, negative 
sequence and zero sequence respectively with taking into 
account the shunt capacitances effect - 
5 As a result of performing the first iteration for the positive 
sequence (Figure 7) the compensated current LM_camp_i is 
calculated and the last index in the subscript denotes the 
first iteration. The calculation is based on deducing the 
shunt current from the calculated positive sequence current l M 
10 calculated from the measured phase currents - Fig. 2: 

L A \_comp_i =Lm -0.5d v lB iL A tasAll V M (12) 
where: 

d v - distance to fault calculated without taking into account 
the shunt capacitance effect (10), 
15 Z - total line length (km) 



tanh 






0.5Z! lL B iL d r l 



B^ L -2£?£uL _ positive sequence admittance (capacitive) of a 
t 

line per km length (S/km) 

Z\ L z=zZ= J ± ' - positive sequence impedance of a line per km length 
20 (£2/km) 



Positive sequence impedance of a faulted line segment, 
between points A and F, without taking into account the shunt 
capacitances effect and considering the lumped line model 
25 equals: 

d v iz\ L (13) 

while for the considered here distributed long line model: 

^/Z'uAinhl {14) 



where : 



Thus, the positive sequence impedance of a line with taking 
5 into account the distributed long line model (Z^ s ) equals: 

Zlong AT (15 ) 

As a result of performing the first iteration for the negative 
sequence {Figure 8) the compensated current I_Ai_comp_\ ls 
10 calculated and the last index in the subscript denotes the 

first iteration. This is based on deducing the shunt current 
from the calculated negative sequence current Lai ■ calculated 
from the measured phase currents - Fig. 2) : 

L A 2_c 0mp _> =La2 -0.5d v lg 2L Au^2Z A 2 (16) 

15 

where, taking into account that the line parameters for the 
positive and for the negative sequences are identical 

4tanh2 = 4tanhl 
20 B 2L =B\ L 

As a result of performing the first iteration for the zero 
J sequence (Figure 9) the compensated current L A o_camp_\ *- s 

calculated, last index in the subscript denotes the first 
25 iteration. This calculation is based on deducing the shunt 
current from the calculated zero sequence current Lao • 
calculated from the measured phase currents - Fig. 2: 
lA0_com P _i = Lao ~ 0-5d v lB^ L A tmh0 V A0 C 17 > 

where : 



tanh 

i^tnnbO ^ 
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O.SZ OL B OL d v l 

»• - J^o/. _ Z ero sequence admittance (capacitive) of a line 
per km length (S/km) 

2^==^- - zero sequence impedance of a line per km length 
5 (Q/km) 

Zero sequence impedance of a faulted line segment, that is 
between points A and F, without taking into account the shunt 
capacitances effect and considering the lumped line model 
10 equals : 

*&l < 18) 

while for the considered here distributed long line model: 

dJZ 0L A s]nh0 <19> 
15 where : 



£±smhQ 



Thus, the zero sequence impedance of a line with taking into 
account the distributed long line model {Z$£ g ) equals: 
20 Z&=A**oZ*l (20) 

The quadratic complex formula (8) with two unknowns 
< d com P 1 tP- u -l - sought fault distance, R F - fault resistance) 
after introducing the compensation (first iteration) takes the 
25 following form: 

<21) 



20 



where : 



10 



15 



20 



ycomp_\ _ —A-P _ f au it loop impedance calculated with: 
La_p 

V_ A p - original fault loop voltage (2}, 

, com P _i _ n r i - fault loop current 

composed of the positive (12), negative (16) and zero (17) 
sequence currents obtained after deducing the respective 
capacitive currents from the original currents. 



Writing (21) in a more compact form results in: 

Ar-'(^..) 2 ^r*- Wi +£ omp - 1 -° 

where : 

at"- 1 = + k:z-' = -mst 1 

^00 _ /1 00_ Re ~*~ J"00_ Im ^ comp _ 1 



(21a) 



z «w«>_i - — a -p fault loop impedance calculated wi 

—A_p jComp^X 
LA_p 

V a - original (uncompensated) fault loop voltage (2) 



with: 



LaZ~ 1 =ZlL A \_comp_\ +a 2 L A z.co m p_i + ^i)LA0_c„ m p_\ 



- fault loop current 
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composed of the positive (12), negative (16) and zero (17) 

sequence currents obtained after deducing the respective 

capacitive currents from the original currents, 

K\. L\> Mi ~ coefficients gathered in TABLE 3 . 

Formula (21a) can be written down separately for real and 

imaginary parts: 
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»W + «- 1+ ^!fc'^= 0 < 21b > 
«Er , W«^-i) a + 4jE" W.i + «"' =o (210 

5 Combining (21b) and (21c) in such the way that fault 
resistance is eliminated, that is, equation (21b) is 
multiplied by and equation (21c) by Afife 1 and then 

subtracting them, yields the quadratic formula for a sought 
fault distance: 

io 2?r p -'(<w,) 2 +Br p - 1 =o < 22 > 

where : 

Br*-* = ^e-X^m 1 -^-'^STr; 1 

15 

Equation (22) has two roots [ {d comp _^ , (d comp l ) 2 1 for the 
distance to fault: 



(4 compel h ?= 



(d cotn p _i ^2 2g com P— 1 



(23) 



■ 20 The root, which corresponds to the selected previously root 

(10) for d (uncompensated) is taken as the valid result. The 
] compensation procedure requires iterative calculations, 

performed until the convergence is achieved (i.e. until the 
] location estimates cease to change from the previous 

* 25 estimates) or as with a fixed number of iterations such as 2 

iterations. The calculated distance to a fault from a 
[ particular (say, present iteration) is utilized for 

determining the shunt current in the next iteration. 



The method of the invention is illustrated in two flow-charts 
of the FL algorithm, Figure 10, single line and in Figure 11, 
parallel lines. 

As shown in the flowchart in Figure 10 the following 
measurements are utilized: 

- voltages from the side A from particular phases a, b, c: 

v A_a ' v A.b ' Ha_c 

- currents from the side A from particular phases a, h, c: i A a 

• *A_b ' 'i4_C 

The input data utilized at step 101, Input data and 
measurements, are as follows: 

- impedances of the line for the positive ( Z, t ) and zero (Zo/,) 
sequences , 

- impedance of the extra link 25, 45, 55 between the 
substations A, B for the positive (negative) sequences (Z lAB ) 

- source impedances for the positive (negative) sequences 

( Z M , Z lsB ) : the representative values or the measured values 
are used, and a communication means is. used for sending the 
measured remote source impedance as previously described, 

- information on the fault type (from the protective relay) . 

The measured fault quantities (voltages and currents) undergo 
adaptive filtering at step 104, Adaptive filtering of phase 
quantities, aimed at rejecting the dc components from currents 
and the transients induced by Capacitive Voltage Transformers 
(CVTs) . 

In the next step the symmetrical components of voltages and 
currents are calculated, step 105, which is equivalent to that 
as shown in Figures 3a, 3b. The fault loop signals are 




composed: the fault loop voltage as in (2), while the fault 
loop current as in (3) - but with taking a 0 m=°- 

The distance to fault without taking into account the shunt 
capacitances effect (d) is calculated at step 106 by solving 
5 the quadratic formula (9). The solution of (9) is presented in 
(10) . 

The result obtained without taking into account the shunt 
capacitances effect d following 106 is treated as the starting 
value for performing the compensation for shunt capacitances. 
10 The distributed long line model is applied for the 
compensat ion . 

The following additional data is required for the calculating 
the compensation for shunt capacitance, step 107: 

- positive sequence capacitance of the line { C XL ) 
15 - zero sequence capacitance of the line ( C 0L ) 

- line length (1) , which is used to express line impedances / 
capacitances per km length. 

The first iteration of the compensation leads to the quadratic 
equation (22), which is solved in (23). The next iterations 
20 are performed analogously. Iterative calculations are 

performed until the convergence is achieved or a fixed number 
of iterations, i.e. 2-3 iterations, may be made. The 
calculated distance to a fault from a particular (say, present 
. . iteration) is utilized for determining the shunt current in 

*25 the next iteration. After completing the iterative 

calculations the distance to fault d camp is obtained. 



As shown in the flowchart in Figure 11 for parallel lines the 
following measurements are utilized: 
V\30 - voltages from the side A and line I»A from particular phases 

- r . 

a, b f c: a , v M b , v M c 



- currents from the side A and line LA from particular phases 

a, b, Ci i/U_a ' *AA_b> l AA_e 

- zero sequence current from the healthy parallel line LBz i^o 



The input data utilized in step 111, Input data and 
measurements, are as follows: 

- impedances of the faulted line for the positive (Z JM ) and 
zero {Zoia > sequences 

- impedance of the healthy line for the positive (negative) 
sequence ( Z ltB ) 

- impedance of the extra link between the substations A, B for 
the positive (negative) sequences 

- zero sequence impedance for the mutual coupling (Z 0m ) 

5 - representative values of the source impedances for the 
positive (negative) sequences 

( Zl\sA ' ) 

- information on the fault type is obtained from the 
protective relay. 

0 The measured fault quantities, the voltages and currents, 
undergo adaptive filtering at step 114 for the purpose of 
rejecting the dc components from currents and the transients 
induced by Capacitive Voltage Transformers (CVTs) . 

in the next step 115 the symmetrical components of voltages 
25 and currents are calculated as shown in Pigures 3a, 3b. The 

fault loop signals are composed: the fault loop voltage as in 
(2), while the fault loop current as in (3) - but with taking 

The distance to fault without taking into account the shunt 
30 capacitances effect (d) is calculated in step 116 by solving 
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the quadratic formula (9) . The solution of (9) is presented in 
(10). 

The result obtained without taking into account the shunt 
capacitances effect (d) is treated as the starting value for 
5 performing the compensation for shunt capacitances. The 

distributed long line model is applied for the compensation. 

The following additional data for the faulted line is required 
for the compensation for shunt capacitances step 117: 

- positive sequence capacitance of the line {C lL ) 
10 - zero sequence capacitance of the line {C 0L ) 

- line length (2) , which is used to express line impedances / 
capacitances per km length. 

In the case of a single line, the compensation is performed 
15 analogously. Iterative calculations are performed until the 
convergence is achieved or by using a fixed number of 
iterations, e.g. 2-3 iterations. The calculated distance to a 
fault from a particular iteration, for example the present 
iteration, is utilized for determining the shunt current in 
20 the next iteration. After completing the iterative 

calculations the distance to fault is obtained. 

Figure 18 shows an embodiment of a device for determining the 
distance from one end, here shown as end A of a section of 
power transmission or distribution line A-B, to a fault F on 
: .-.25 the power line according to the method of the invention 

: described. The fault locator device 1 receives measurements 
..!.: from measuring devices located at one end A such as current 
,.J measuring means 14, and voltage measurements from voltage 
.-;-*. measurement means 11. The fault locator device may comprise 
/...:30 measurement value converters, members for treatment of the 

calculating algorithms of the method, indicating means for the 
calculated distance to fault and a printer or connection to a 
.-•-! printer or facsimile machine or similar for printout of the 
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calculated fault- In a preferred embodiment of the device the 
fault locator comprises computer program means for providing a 
display of the information provided by the method of the 
invention, such as distance to a fault d or d comp on a terminal 
which may be remote from the location of the line and/or the 
fault locator. Preferably the computer program means receives 
information from the fault locator and makes it available to 
provide information for a display of a computer such that an 
operator or engineer may see a value for the calculated 
distance to a fault displayed- The value may be displayed 
relative to a schematic display of the line or network in 
which the fault has taken place. 

In the embodiment shown, measuring device 14 for continuous 
measurement of all the phase currents, and measuring device 11 
for measurement of voltages are arranged at one end, station 
A. Optionally, measuring devices such as 15, 13 may also be 
arranged at station B but they are not necessary to practise 
the invention. The measured values such as: the three-phase 
voltages (VXa) of the faulted line, three-phase voltages ( Iaa) 
of the faulted line and zero sequence current {l AB o) from the 
healthy parallel line (note that zero sequence is not present 
when the single line only is considered) , and a value 
representative of the source impedance at B, Z UB as a are all 
passed to a calculating unit comprised in fault locator 1, 
filtered such as described in relation to figure 3a, 3b, and 
stored in memory means. The calculating unit is provided with 
the calculating algorithms described, and programmed for the 
processes needed for calculating the distance to fault. 
Optionally, the source impedance for the remote end, Z UB may 
be measured by the remote device RD, 10, and the information 
sent via a high speed communication means 60 to the fault 
locator at A. In some applications it will be preferable to 
use a measured value sent from B instead of a representative 
value stored at A. It may be seen in Figure 18 that current 



measuring means 15 and voltage measuring means 13 at remote 
end B may provide a RD 10, a fault locator, or any suitable 
device with measurements to calculate the remote source 
impedance 

The calculating unit of fault locator 1 is provided with pre- 
fault phase currents and also with known values such as shunt 
capacitances and the impedances of the line. In respect of the 
occurrence of a fault, information regarding the type of 
fault, phase-to-phase, phase- to- ground etc., maybe supplied 
to the calculating unit of the fault locator. When the 
calculating unit has determined the distance to fault, it is 
displayed on the device and/ or sent to display means which may 
be remotely located. A printout or fax of the result may also 
be provided. In addition to signalling the fault distance, the 
device can produce reports, in which are recorded measured 
values of the currents of both lines, voltages, type of fault 
and other measured and/or calculated information associated 
with a given fault at a distance. 

The method and a fault locator device according to any 
embodiment of the invention may be used to determine distance 
to a fault on a section of power transmission line. The 
present invention may also be used to determine a distance to 
a fault on a section of a power distribution line, or any 
other line or bus arranged for any of generation, 
transmission, distribution, control or consumption of 
electrical power. 

The fault locator device and system may comprise filters for 
filtering the signals, converters for sampling the signals and 
one or more micro computers. The microprocessor (or 
processors) comprises a central processing unit CPU performing 
the steps of the method according to the invention. This is 
performed with the aid of a computer program, which is stored 




in the program memory. It is to be under stood that the 
computer program may also be run on one or more general 
purpose industrial computers or microprocessors instead of a 
specially adapted computer. 

5 

The computer program comprises computer program code elements 
or software code portions that make the computer perform the 
method using equations, algorithms, data and calculations 
previously described- A part of the program may be stored in a 

10 processor as above, but also in a ROM, RAM, PROM or EPROM chip 
or similar. The program in part or in whole may also be stored 
on, or in, other suitable, computer readable medium such as a 
magnetic disk, CD-ROM or DVD disk, hard disk, magneto-optical 
memory storage means, in volatile memory, in flash memory, as 

15 firmware, or stored on a data server. 

A computer program according to the invention may be stored at 
least in part on different mediums that are computer readable. 
Archive copies may be stored on standard magnetic disks, hard 
20 drives, CD or DVD disks, or magnetic tape. The databases and 

libraries are stored preferably on one or more local or remote 
data servers, but the computer program and/or computer program 
product may, for example at different times, be stored in any 
of; a volatile Random Access memory (RAM) of a computer or 
25 processor, a hard drive, an optical or magneto-optical drive, 
or in a type of non-volatile memory such as a ROM, PROM, or 
EPROM device. The computer program may also be arranged in 
part as a distributed application capable of running on 

: several different computers or computer systems at more or 

: 30 less the same time. 



In another preferred embodiment, the fault locator may be used 
with parallel lines to locate single phase- to-ground faults 
<a-g, b-g, c-g faults) in case of unavailability of 



measurements from the healthy parallel line. Two modes of the 
healthy line operation are taken into account: 

- the healthy line being in operation, 

- the healthy line switched-off and grounded at both the ends. 

Figure 19 shows the flow-chart of the algorithm for locating 
faults in parallel lines under unavailability of measurements 
from the healthy parallel line- The unavailable, here healthy 
line zero sequence current, is required for reflecting the 
mutual coupling effect under single phase-to-ground faults (a- 
g, b-g, c-g faults) , The unavailable current is thus 
estimated. The other faults can be located with the standard 
fault location algorithm (such as the algorithm from reference 
[1]) • 

The sequence of computations for the presented one-end fault 
location algorithm is as follows. 

As shown in the flowchart in Figure 19 the following 
measurement s are utilized: 

- voltages from the side A and line LA (faulted) from 
particular phases: , v /U _ & , v^_ c 

- currents from the side A and line LA (faulted) from 
particular phases: Iaa^o » *aa_*' *aa_c 

The utilized input data are as follows: 

- impedances of the faulted line for the positive {Z llA ) and 
zero (Z 0LA ) sequences 

- impedance of the healthy line for the zero sequence {Z QLB ) 

- zero sequence impedance for the mutual coupling <Z 0m ) 

- information on the fault type (from the protective relay) 

- mode of the healthy line operation: in operation or 
switched-of f and grounded at both ends) - 

The measured fault quantities (voltages and currents) undergo 
adaptive filtering aimed at rejecting the dc components from 



currents and the transients induced by Capacitive Voltage 
Transformers (CVTs) . 

The following equations are used in the method of the present 
parallel line embodiment. In addition to the algorithm (1) 
described above for estimating distance to fault (d [pu] ) by 
considering the Kirchhoff's voltage law for the fault loop as 
seen from the locator installation point: 



R ( Zfi^Iaai x ZfiLaai . v *fqLaa*\_q (24) 
^ Kfi %-F2 £fo ) 



the fault loop voltage (VaajJ and current (Laa_j>) can be 
expressed in terms of the symmetrical quantities: 

and 

/ _ r SL i ^Qm / (25) 
where : 



iJh .=^z..,+^/^+flo#^/>uo < 25 *> 

\LA 



LAA p = + *2±./U2 + "0 *^ ±AAQ 



is the fault, loop current without* compensating . for the mutual 
coupling effect (i.e. composed as for the single line - 
superscript SL) , 

«2' £o~ complex coefficients gathered in TABLE 1 (the 
derivation as in APPENDIX Al) , 

Y-AAi' Y.AA2' Y.AA0 ~ positive, negative and zero sequence of 
measured voltages, 

Laa\* Laa2' Laao - positive, negative and zero sequence current 
from the faulted line LA, 

Labq - unavailable zero sequence current from the healthy 
parallel line LB (to be estimated), 




ZiiA' £qia ~ positive and zero sequence impedance the whole 
line LA, 

Z Qm - zero sequence impedance for mutual coupling between the 
lines LA and LB, 
5 R F - -unknown fault resistance. 

In the next step the symmetrical components of voltages and 
currents are calculated as shown in Figures 3a, 3b, The fault 
loop signals are composed: the fault loop voltage as in (2) , 
while the fault loop current as in (25). The formulae (2)-(25) 
10 present the fault loop signals expressed in terms of the 

symmetrical components of the measured signals. However, one 
can use the classic way for composing the fault loop signals . 



The presented method covers single phase- to ground faults (a- 
15 g, b-g, c-g faults) . The other remaining faults have to be 

located with the fault location algorithms described above or 
standard fault location algorithm [1] . Distance to fault (d) 
for the considered here single phase-to ground faults is 
calculated by solving the quadratic formula for a sought 
20 distance to .f ault . (26) • Equation (.2.6) is the .same as. equation 
(10) except that the values for B l9 B 2 ,B 3 are different to the 
values determined in (10). The solution gives two roots: 

2B * (26) 

d ^ 25 

(as above, the root which fulfils the condition (0<d<l) is 
25 selected as the solution for the distance to fault) . One has 
to substitute the following into (26) : 
B 2 - real(A 2 ) imag(A 00 ) - realC^ ) imagCA 2 ) 
B { =real(^,)imag(A 00 )--real(A 0(> )imagCA 1 ) 



B 0 ==real(A 0 )imag(A 00 )-real(A 00 )iniag(A 0 ) 
where : 

Ai =~Z\[aKiLaa_p — —ILiLaao - p ^ i£ 0 

£-0 

The recommended SET of the coefficients b F1 ,b F2 are taken from 
10 TABLE 4 and the recommended SET of the coefficients a Fl , a F2 , 
a F0 from TABLE 5. 

Fault loop voltage in this embodiment is found from the TABLE 



below 



Fault loop voltage composed 
in terms of symmetrical 
components 


Fault loop voltage composed 
as in the classic approach 


V AA ^ p =a x V_ AM + a 2 V AA2 +a 0 V MQ 
a-g fault: a x ^a 2 —a 0 =l 
jb-g fault : flj = fl 2 , a 2 —a, a Q = 1 
c-g fault: a x =a, a 2 = a , a Q = 1 
a = exp(j2#/3) , j = V =: T 


a-g fault: V M _ p =V^_ fl 
a-g fault: K/u_ P =Ka/\_* 
a-g fault: Za*_p = £/w_c 



15 

Fault loop current I^_ p composed as for the single line is 
found from the TABLE below 



# 



33 



rSL 



in terms of symmetrical 



components 



Z 

LaA p = £\Iaa\ + «2lA42 +50 ~ LaAQ 

a-g fault: a j = £ 2 = = ^ 

Jb-fir fault: =£ 2 , a 2 =a / fi 0 =1 

c-gr fault: £1 / £2 = ^ » £0 = ^ 
a = exp(j2^r/3) , j = 4-i 



as in the classic approach 



a-gr fault : Laa_p = Ia4_« + £o£a40 
a -fir fault: Zj£_ p +£o£mo 
a-gr fault: I^_ p = I>w_<2 ^^oLaaq 
where: fc 0= ^=^ 
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The complex coefficients dependent on the mode of the healthy 
parallel line operation: 

a) healthy line LB is in operation: 



Ho 



M , =Z 1M Z ILfl + Z lw (Z IxA +Z hpB ) + Z l£jB (Z UA +Z, yB ) 



15 



b) healthy line I/B is switched-off and grounded: 




in another embodiment, which is presented here, the method is 
applied for the standard availability of the measured signals 
and is only valid for ground faults, including both: 

- single phase-to-ground faults 

5 - phase-to-phase-to-ground faults. 

Fault location procedures obtained for these faults under the 
standard availability of the fault locator inpufc signals are 
extremely simple and compact. Distance to a fault is 
calculated with a first order formula. 

10 Figure 22 shows a flow-chart of the developed algorithm for 
locating ground faults in parallel transmission lines. 
The sequence of computations for the presented one-end fault 
locator is as follows. As shown in the flowchart in Figure 
22. B3 the following measurements are utilized: 

15 - voltages from the side A and line LA from particular phases 

a, br c: v M _ a t v M _b » v M c 

- currents from the side A and line LA from particular phases 

a, by C: lAA_ a , *AA_b ' ^AA_c 

- zero sequence current from the healthy parallel line LB: i AB0 

20 The utilized input data are as follows: 

- impedances of the faulted line for the positive (Z 1LA ) and 

zero ( Z QlA ) sequences , 

- impedance of the healthy line for the zero sequence (Z 0LB ) 
; - zero sequence impedance for the mutual coupling (Z 0m ) 

25 - information on the fault type from the protective relay. 

: The measured fault quantities (voltages and currents) undergo 

: adaptive filtering aimed at rejecting the dc components from 

: currents and the transients induced by Capacitive Voltage 

: Transformers (CVTs) , preferably as described and shown in 

:30 relation to Figures 3a, 3b. 
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A generalized model of the fault loop, used for deriving the 
algorithm of the present embodiment, is stated as follows: 

V M _ p - dZ llA l_AA_p ~ R F i&FiLFX + ZfiLf 2 + 3. fUfo) = 0 ( 1 > 

where : 

5 d - unknown and sought distance to fault, 

Z-iiA ~ Positive sequence impedance of the faulted line, 
V Mp ,Laa p - fault loop voltage and current composed according 
to the fault type, 
R F - fault resistance, 
10 Lr, ~ sequence components of the total fault current (i=0 - 
zero sequence, i=l positive sequence, 

i=2 - negative sequence) , 
a Fi - weighting coefficients (TABLE 2) . 

Fault loop voltage and current can be expressed as in classic 
15 distance protection technique or, as in this document, in 
terms of the local measurements and with using the 
coefficients (£ 0 < «i . &2> which are gathered in TABLE 1 
(derivation of the coefficients is contained in APPENDIX 
APP.l) : 

L M _ P = *LMl + ^L M 2 + «0 + -°^ LaB ° 1 3 ) 

where : 

AA, AB - subscripts used for indicating measurements acquired 
from the faulted line (AA) and from the healthy line (AB) , 
25 respectively, 
: Zom ' Z-om " impedance of the faulted line and mutual coupling 

between the lines for the zero sequence, respectively. 

In the next step the symmetrical components of voltages and 
'*: currents are calculated as shown in Figure 3a, 3b. The fault 



loop signals are composed: the fault loop voltage as in (2), 
while the fault loop current as in (3). The formulae (2) -(3) 
present the fault loop signals expressed in terms of the 
symmetrical components of the measured signals. However, one 
5 may instead use the classic way for composing the fault loop 
signals, as shown in APPENDIX Al. 

The presented method covers single phase- to ground faults (a- 
g, b-g, c-g faults) and phase- to -phase -to-ground faults (a-b- 
g , fc>-c-g, c-a-g faults) , thus, the faults for which the 
10 highest fault resistance can be expected. The other remaining 
faults have to be located with the fault location algorithms 
described above or a standard fault location algorithm, such 
as for example the fault locator from reference [1] ) . 

Distance to fault (d) for single phase- to-ground faults is 
15 calculated as follows: 

a-g fault: 

imag{VAA^(lAA0 -PqLabq)F} (27a) 



where : 



Fault loop signals composed 
in terms of symmetrical 
components 


Fault loop signals composed 
as . in the classic approach 






lAA_p -&\LaA\ +£21aA2 + 

+ £o - Laaq 1 2o - Iabo 

£l\LA 4-lLA 


LAA_p ~IaA_€i +kolAAO+!tOmlABO 


where : a x = a 2 = a Q = 1 


where : k 0 - — - . , fc Qm - 

AiLA *±\LA 



P Q = ~ 0LB ~ 0m (for symmetrical lines: P 0 = 1 ) . 




b-g fault: 



• 


j^^^^. t * 

■ H : - • • * 

• — - . : 


imag{V M B [3a 2 aMO-£olA5o)]*) 


(27b) 


*mag { (Z IM J aa _ p )[3^ a aao ~ EoLabo )r 1 
where : 


Fault loop signals composed 
in terms of symmetrical 
COluponento 


Fault loop signals composed 
as in the classic approach 






+So „ Iaao + So ~ — Iabo 


Lm„P =!-AA_b +!$-oLaAO + *OmlAB0 


2 _ 1 
wuere : o j ££ » ££2 ££ * — o * 

a = exp(j27T/3) 


wnere : «o" ~. * ^0m ^ 


p 0 = — 0/J? = 0m < for symmetrical lines : £ 0 = 1 > • 

^OW ~2om 


c-g fault: 

zmag{Vl M p [3^(7^40 -PqIabo)!*) 


(27c) 


" /ma^{(Z I£-4 /^_ p )[3fl(/ yW0 -PoIaboW) 


where :• 




Fault loop signals composed 
in terms of symmetrical 
components 


Fault loop signals composed 
as in the classic approach 


V M _p =g. { V AAi + £ 2 K/W2+SoKaao 




2oL>* t i _ 2.0m f 

+ So ~ I/wo + So - Labq 


LaA_p -lAA_c +KqLaA0 +k.0mlAB0 


2 ! 

where : a } = a , a 2 = £ ' So 3=5 1 
a = exp(j2?r/3) 


ZlOLA —Z.UA i, — 2.0m 

where: k Q - " ■ * *0m 7 

Zl\LA ^im 



P 0 = - 0LB - 0m (for symmetrical lines: P 0 =1 >- 
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Distance to fault (d) for phase-to- phase-to ground faults can 
be calculated in two different ways, depending whether the 
pre-fault currents can be used or have to be avoided. 
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1 . Procedure for distance to fault calculation with the use of 
pre-fault measurements: 

a-b-g fault: 

imag{(Z llA LAA_ P M.(lAAO-£oLABo)T} 



(28a) 



where: W =MM1±ML6A1 , 
where, w cAL^+W-aaz 



A = l-£ , 



B = l-a , C = l + a 2 , D = l + a 



Fault loop signals composed 
in terms of symmetrical 
c omponent s 


Fault loop signals composed 
as in the classic approach 




V M _ p = V M(l ~^AA„b 


LaA_ P ="iIaM +Q.21AA2 + 

Z.GLA j , ^ £-0m r 
+ LAAO+aO 7 I ABO 

Z\la Mla 


l-AA_v -LAA^a ~I-AA_b 


where : a x = 1 - a 2 , a 2 = 1 - a , d 0 -0 
fi = exp(j2^/3) 





P Q = ~ 0LB — ^Ssl (for symmetrical lines: £ 0 = I > 



15 Jb~ c-gr fault: 

imagiVAA pWiLAAO -EoLabo)Y) (28b) 
~ imag[{Z lLA LAA^MlAAO -EoLabo)Y) 
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, ~. AALaai+MLaai 

where: W=-— —=- 

QALkA\ +£LLaa2 



A = a 2 -a , B = a-q 2 , C = -l , D--1 



Fault loop signals composed 
xn terms of symmetrical 
components 


Fault loop signals composed 
ae Art t*Vif» rl acjcjir aryorosch 


V M? =£,V^i + £2^2 +9*Y-aao 


HaA_P -^AA.fr ~KaA_c 


Im.p = ZiLaai +£iLaai + 

+ *0 ~ Z/UO+SO ~ J- ABO 
Z.UA ^\IA 




where: a x =a 2 -a, a 2 =a-a 2 , fi 0 =0 
a = exp(j2^r/3) 





P 0 = — 0LB - — ^Ssl (for symmetrical lines: £ 0 = 1 > 



c-a-g fault : 

d= imag {V_aa _ p fed, aaq " PqLabo)T ) 
where : 



w _ AALam +M.AA1 
A = a-1, B=« 2 -l, C = a + 1, £ = a 2 +l 



Fault loop signals composed 
in terms of symmetrical 
components 


Fault loop signals composed 
as in the classic approach 


YaA_ P =SjK>UI + *2^AA2 +S 0 ^AAO 




IaA_p = £i!aA1 +£2^AA2 + 

+ *0 - lAAO+fiO ~ —ABO 
£HLA — 1LA 


Iaa_p ~Laa_c ~Iaa„o 


where: ^ = a -1 , a 2 = a 2 ~ * * So = 0 
S = exp027T/3) 





\ : ;- : -~- - ---- 

*■ — -i - - t . 

p 0 = — 0LB 0m (for symmetrical lines: £o=D. 

2 . Procedure for distance to fault calculation without the use 
of pre- fault measurements: 
5 a-b-g fault: 

imaglQLn +V*)(I»ao -£oIabo)*] (29a) 

" imag[Z llA (± a +h+ UoLaao + 2£ 0m I*flo)(I/uo " £o£wo)*l 
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h-c-g fault: 

imagKVfc +H c )(Ii»in -EoLabqW (29b) 

~ imag[Z lLA (± b +L C + 2 KqLmo + 2fcomi-ABo)(I/UO -£oIabo)*1 

c-a-gr fault : 

imag[(V c +V„)(/>un -EoLabq)*} (29c) 

~ imag[Z lIA (.lc + L + 2 !LoLaao + 2fc 0 m2.ABoX£>uo -£oIabo)*1 



where: fc, ^"^ 1 " 

15 




It is also noted that while the above describes exemplifying 
embodiments of the invention, there are several variations a: 
: 20 modifications which may be made to the disclosed solution 
": * without departing from the scope of the present invention as 
defined in the appended claims. 
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Tables 

TABLE 1. Coefficients for determining the fault loop voltage 

p> current U^p) in terms of symmetrical components 

as defined in (2) and (3). 



Fault type 


£i 


£2 


£0 


a-g 


1 


1 


1 




a 2 


a 


1 


c-g 


a 


a 2 


1 


a-b, a-b-g 
a-b-c, a-b-c~g 




\-a 


0 


b-c, b-c-g 


a* -a 




0 


c-a, c-a-g 


a -I 


a 2 -\ 


0 


a = expG2#/3), j = V-l 



TABLE 1A. Fault loop voltage (V >w _p) and current (Laa_ v ) 
expressed with using the classic approach. 



Fault type 


Yaa_p 




a-g 






b-g 




lAA_b + KoLaAQ Oml ABO 


c-g 


YaA„c 


LaA_c + KqLaAO + KomlABQ 


a-Jb, a-b~g 
a-b~c, a-b~ 
c-g 


Y-AA_a ~Y-AA_b 


LAA_a ~~lAA_b 


b-c, b-c-g 


Y^b-Y-AA^c 


I-AA^b ~IaA_c 


c-a, c-a-g 


YaA^c ~YaA_o 


LaA_c ~ l~AA_a 


a = exp02?r/3) # j = V = T 



10 TABLE 2. Alternative sets of the weighting coefficients for 
determining a voltage drop across the fault path resistance 
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Fault 


Set I 
( recommended ) 


Set XI 


Set III 




type 


S.F1 


" f 2 




CL r*t 


— •/* z 




a F1 


&F2 




a-gr 


0 


3 


0 


3 


0 


0 


1,5 


1,5 


0 


b-g 


0 


3a 


0 




0 


0 


1,5a 2 


1,5a 


0 


<- id 


0 


3a 2 


0 


3a 


0 


0 


1,5a 


1,5a 2 


0 


a-Jb 


0 


l-a 


0 




0 


0 


T rr 

0,5(l -a 2 ) 


0,5(1 -a) 


0 


Jb-c 


0 


2 

a-a 


U 


2 „ 




V 


H <Li 2 _ ft \ 


n si/? — /i 2 1 


0 


c-a 


0 


a 2 -I 


0 


a-l 


0 


0 




0,5^ 2 - 1) 


0 


a-b-gr 


I — a 2 


l-a 


0 


l-a 2 


l-a 


0 


l-a 2 


l-a 


0 


b-c-gr 


a 2 -a 


a-a 2 


0 




a-a 2 


0 


a 2 -a 


a-a 2 


0 


c-a-gr 


a-l 


a 2 -l 


0 


a-l 


a 2 -1 


0 


a-l 


s 2 -i 


0 


a-b-c-g 
(a-b-c) 


I- a 1 


0 


0 


l-a 2 


0 


0 


l-a 2 


0 


0 




a = exp(j2;r / 3) , j = -\PT 



TABLE 3. Coefficients for determining a fault current 
distribution factor (6) 



SINGLE LINE (Fig.4) 




Ki - -Z.llZ.UB~iS.lsA +Zi,b)Z\L 

L\ ~Z\lSZ\)A +Zub)+Z\AbIZ.VL + Zub) 

Ml =(Zl»l + Zub)&1AB +ZlLhZ lL ZiAB 


Z.XAB 


K_\ ~~Z\i L\ -Zil + ZisB 

M.1 =Z\tA +Zl,B + ZlL 



PARALLEL LINES (Fig.5) 

K\ =-z, M (z, jA +z lffl +z 1JLB&A5 ) 

Li = z 1LA (z ljA ^Z^g +z 1LB& ^)+z lLS&AB z ljrB1 

Mi =Zaia&\ib*.ab +z iIA (z lsA + z UB )+z {IM&AB (z lsA + z, Jfl ) 

, ~ ?L\LB?L\AB 

where: Z lLBStAB « -T — 




Table 4 The recommended SET of the coefficients b Fl , b F2 i 
relation to (26) 



Fault 




£f2 


a-g 


0 


1 


b-g 


0 


fl 2 


c-g 


0 


a 


a = exp02?r/3) , j = 





Table 5 The recommended SET of the coefficients a FX , a F2 , a F0 
5 in relation to (26) 



FAULT 




«F2 


&F0 


a-gr 


0 


3 


0 


b-g 


0 


3a 


0 


c-g 


0 


3a 2 


0 


a = exp(j2fl7 3) , j = *J-l 
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APPENDICES 

Al. DERIVATION OF THE COEFFICIENTS a,, a 2 , a 0 (TABLE 1) 
Single phase- to-ground fault: a-g fault 

Vaa_ p =V Ma =V AAi +Vaa2 +Y.M0 = <LlY-AM +52^2+^0^0 
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Im_p =LAA_a+hLAAO + k 0m L AB o = Lam + Laai +Lmo + Z0lA z^ llA Laao + zHZ Um " 



Thus : £1 =£2 = ^0 = * 



5 inter-phase faults: a-b, a-b-gr, a-b-c, a-b-c-g faults 
= (l-a 2 + (1 - o)!Laai = 2l £ vui + S2V.AA2 + SoE/MO 

Thus: a x =l-a 2 , a 2 = \-a, a 0 =0 

10 A2. Derivation of coefficients fl Ft > 4f2 > Zfo (TABLE 2) 

TABLE 2 contains three alternative sets (Set I, Set II, Set 
III) of the weighting coefficients, which are used for 
determining a voltage drop across a fault path. The 
coefficients are calculated from the boundary conditions - 
15 relevant for a particular fault type. This is distinctive that 
in. all the sets the zero sequence is omitted ( a F0 =Q) . This is 
advantage since the zero sequence impedance of a line is 
considered as the uncertain parameter. By setting a F0 =0 we 
limit adverse influence of the uncertainty with respect to the 
20 zero sequence impedance data upon the fault location accuracy. 
To be precise one has to note that this limitation is of 
course partial since it is related only to determining the 
voltage drop across a fault path. In contrast, while 
determining the voltage drop across a faulted line segment the 
25 zero sequence impedance of the line is used. 



a-gr fault, Figure 12: 

30 



Taking into account that in the healthy phases: l F b - l F c - 0 
what gives: - I_ F , =|fc F _ fl +«If_6 + <? L F _c)=\kp _a + <*>+° 2 «)=\Lf_* 

The sequence components are related: Lpi =L F2 =L F0 311(1 
finally : 

I F =l F _ a =3L F 2> thus: fifi=°» fiF2=3, fifo=° * as in the Set x . 
from Table 2) 
or 

J F = / F _ fl = 37 F1 , thus: fi F1 =3, S F2 =0, a F0 =0 (as in the Set II 
from Table 2) 



or 



7 F =I f _ fl =l,5/ FI +l,5/„, thus: fiK1 =l,5, a F2 =U, a F0 =0 (as in the 
Set III from Table 2) 

a-jb fault Figure 13a, 13b: 

The fault current can be expressed as: L F =L F _ a or as : 

Taking into account that in the healthy phase: Lf_c=° and 
for the faulted phases: L F _b=-lF_o> what gives: 

LFO = \k F _« +L F _ b +hJ = \ll F _<,+{-lFj+<>h Q 

The relation between l FI and l F2 is thus: 
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Finally: 

If =L F _ a = -^Z^) Lf2 
thus: a f7 =0, a F2 £fo=<> < as in the Set 1 from Table 2) 



or 



thus: a F1 =l-o 2 , S«=0- Sfo=0 ^s in the SET 2 from Table 2) 
or 

£ f = 0,5/ , _ a + 0.5/ F _ fl = -^pjLn + ^Ifi - 0.5(1 - a)/ , 2 + 0,s(l - a 2 )l fi 

thus: a Fl = 0,5(l -« 2 ), a F2 =0,5(l-a) , a F0 =0 (as in the Set III 
10 from Table 2) 

(a-b-g) fault. Figure 14: 

* Lf = Lf_« ~lF_b = (Lf\ +Lfz +Ifo)-(« 2 Ifi +«If2 +Lfo)= 

15 Thus: a Fl =l-a 2 , a F2 = l-a, a F0 = 0 (as in the Sets I, II. Ill 
from Table 2) 

(a-b-c) or (a-b-c-g> symmetrical faults. Figure 15a, 15b, 15c: 

20 Taking the first two phases (a, b) for composing the voltage 
drop across a fault path one obtains: 

Lf =Lf_o -L F _b = Ufi +Lfz +Lfo)-((L Z Lfi +«If2 +Lfo)= 

= (l-a 2 )/ n +(l-£)l F2 
Thus : 

a pi = 1 - a 2 , a F2 = 1 — a , a F0 -0 
25 Additionally, if a fault is ideally symmetrical the positive 
sequence is the only component, which is present in the 
signals. Therefore, we have: 



a Fl =l-a 2 , a F2 =0, a F0 =O (as in the Sets I, II. Ill from Table 
2) . 

A3. Derivation of the complex coefficients in the fault 
current distribution factors for the positive (negative) 
sequence (Table 3) 



a) Case of the single line with extra link 45 between the 
substations (Figure 4) 

Let us determine the fault current distribution factor for the 
positive sequence (the fault current distribution factor for 
the negative sequence is the same) . The equivalent circuit of 
Figure 4 with indicated flow of currents for- the incremental 
positive sequence is presented in Figure 16. 

Considering the closed mesh containing: the local segment of 
the faulted line, the remote segment of the faulted line and 
the extra link between the substations, one can write down: 

dZ lL AL M + (1 - d)Z lL (AL M -LfO-ZiabALabi = & 

From the above equation the unknown current from the extra 
link between the substation can be determined as: 

Considering the closed mesh containing the source impedances 

£t*fl> 3X143 the extra link ( - ,Afl) one Can write down: 

Zua(4Ui +^LabO + Z 1ab AI_ abi +Z UB (AL Ai +AL ABl -Z.fi) = ° 

Substituting the previously determined unknown current from 
the extra link into the above equation one obtains: 
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where, as in TABLE 3 (Single line, Z MB *~ ) , we have: 
5 L 1 =Z il (Z m + Z 1 ,b)+Z, ab {Z il + Z 1sB ) 

If there is no extra link between the substations (Z MB ->«) 
has to consider the closed mesh containing the source 
impedances (Z isA , Z lsB ) and both the segments of the faulted 
10 line ldZ XL and Q-d)Z lL 1 . For this mesh one can write: 

[Z XsA +dZ lL )AL Al +\Z UB +{\-d)Z lL ](AL M -Ifi) = ° 
After rearrangements one obtains: 

where, as in TABLE 3 (Single line, Z iAB -*« ) , we have: 
15 £, =~Z lL 

Mi = z UA +z isA + z lL 



b) Case of the parallel lines with extra link between the 
20 substation (Figure 5) 



Let us determine the fault current distribution factor for the 
positive sequence (the fault current distribution factor for 
the negative sequence is the same) . The equivalent circuit of 
25 parallel lines from Fig. 5 with indicated flow of currents for 
the incremental positive sequence is presented in Figure 17. 
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The healthy parallel line [LB) and the extra link 55 (AB) , 
which are in parallel connection, have been substituted by the 
equivalent branch with the equivalent impedance: 

£\IB + ^-lAB 

5 considering the closed mesh denoted by (AA, F, BA, BB, AB) one 
can write: 

From the above equation the unknown current from the 
equivalent branch can be determined as: 

10 ALmabi =^ IiM -(4Lmi -V-VL F1 ) 

B.XLB&.AB 

Considering the closed mesh containing the source impedances 
{Z XsA , Z XsB ) and the equivalent branch (Z il£&AB ) one can write 
down: 

Z\sa(ALaa\ +ALlb&abi) + z -\lb&abALlbslabi +Z.i*b(. a Lam + a Llb&abi -In) = ° 

15 Substituting the previously determined unknown current from 
..the healthy line- into the above equation one obtains : 



Lfx 



where, as in TABLE 3 (Parallel lines), we have: 
K\ = ~Z XLA {Z UA + Z Jjfl +Zl\lb&.ab) > 

Mi = z xlA z lLB&AB + z ilA (z XsA +z XsB )+z lLBScAB (z UA +Z UB ) 

, ~ Z\lbZl\AB 

where : A\lb&ab - ~ T~y 

In case if the extra link between the substations ( Z XAB ) is 
not present one has to substitute: Z XLB&AB =Z 1LB . 

25 
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CLAIMS 

1. A method to locate a fault from one end of a section of a 
power line (A-B) by means" of measurements of current , voltage 
and angles between the phases at a first (A) end of said 
5 section, and that upon detection of a fault condition between 
said first end and a second end of said power line, 
characterised by 

-calculating symmetrical components of currents for said 
current and voltage measure at said first end, 
10 -calculating a distance (d) from said first end (2) to the 
fault (F) the distance (d) to the fault using a quadratic 
equation of the form: 
B 2 d 2 +Bid + B 0 =0 
where : 

15 B 2 ~ A 2 _ReA>0_Im _ ^2_ImA)0_Re 

B x = A iR( .A m i m — A 1 _i m i4 o0Re 

B 0 = A)_ReA>0_Im ~ 4>_Im4)0_Re 

2. A method according to claim 1, characterised by calculating 
20 the distance (d) to the fault using an equation of the form: 

—AA _Jp 

where : 

z = — AA -?- - calculated fault loop impedance. 

25 3. A method according to any of claims 1 or 2 , characterised 

by calculating the distance (d) to the fault using an equation 

of the form: 

A 2 d 2 + A x d + A 0 + A 00 jR f " 0 

where : 

30 A 2 = A 2 Re 4- jA 2 _ Im = K x Z llA 
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A 0 = i4o_ Re + jA)_im - -Li2Laaj 

Al0_ Re + -K*O0_Im — />u_p 

_ Y-AAjp _ calcu i a ted fault loop impedance 
Laa_p 

Ki. U> M.i = coefficients gathered in TABLE 3 



4. A method according to any of claims 1-3. characterised by 
-determining source impedance at said first end as a 
representative value, and 

-determining a value for source impedance at said second end 
10 as a representative value. 

5. A method according to any of claims 1-4, characterised by 
calculating a value for impedance of an extra link (45, 55) 
between the terminals A, B, as having impedance for the 

15 positive sequence equal to: 
, 7 Z lLB Z [AB 

I 4l\LBSlAB ~ - ? > 

<L\LB + ^MS 

where 

Z lAB = impedance for the positive sequence of the extra link, 
Z 1M = positive sequence impedance of the healthy line. 
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6. A method according to any of claims 1-5, characterised by 
calculating symmetrical components of currents for said 

' : current and voltage measured at said first end by: 

-inputting instantaneous phase voltages (30a) , 
25 -filtering (33a) the values to determine the phasors, and 

-: -calculating (34a) phasors of symmetrical components of 

voltages . 

7. A method according to any of claims 1-6, characterised by 
A 30 calculating symmetrical components of currents for said 

*: current and voltage measured at said first end by 
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-inputting instantaneous phase currents and instantaneous zero 
sequence current from a healthy line (30b) , 
-filtering (33b) the values to determine the phasors, and 
-calculating (34b) phasors of symmetrical components of 
currents . 

8. A method according to any of claims 1-7, characterised by 
determining a compensation for shunt capacitance by means of 
an equation of the form: 

jt^gw.!? +Br mp - l d comp _ l +B^- i =o <22) 

where : 

Br"- 1 = A c :^-%Tta ! - a"z-*ast& 

9. A method according to claim 8, characterised by 

determining a compensation for shunt capacitance by means of 
an equation of the form: 

ar'-'^w.) 2 +4r p - , <wi + ^ mp - i =o (2ia> 

where : 

7 compel _ ¥-A-p fault loop impedance calculated from: 

V a - original (uncompensated) fault loop voltage, 

/rr , =^i-c^_i + fi2lA2_ co „ v _»+aoI>io_ e ^_. - ^ult loop current 
composed of the positive (12), negative (16) and zero (17) 




sequence currents obtained after deducing the respective 
capacitive currents from the original currents, and 
K\> Ll\* Mi = coefficients gathered in TABLE 3. 

5 10. A method according to any of claims 1-9, characterised by 
measuring the source impedance Z l)A at said first end A. 

11. A method according to any of claims 1-9, characterised by 
-measuring the source impedance Z UB at said second end B, 

10 -sending a communication of the measured value of source 

impedance Z isB at said second end B to a fault locator at said 
first end A. 

12. A method according to any of claims 1-11, characterised by 
15 -determining the zero sequence current from the healthy line 

of a section of parallel power lines, 

-calculating a distance to a fault for the parallel line 
section. 

20 13. A method according to claim 12, characterised by 

determining distance to a single phase to ground fault without 
measurements from an operating healthy parallel line by means 
of complex coefficients P 0 according to a formula of the form: 

2 5 and K x , L, , M, according to 

Mi ~?l\la^\lb + z IM (z lx4 +z UB )+z iLB (z hA +z isB ) 

3 0 14. A method according to claim 12, characterised by 

determining distance to a single phase to ground fault without 




measurements from switched off and grounded parallel line by 
means of complex coefficients P 0 according to 

and K lt L, , M_ , according to 
5 K\ — ~Zl\la 

Ll = Z.ILA + £lsB 



15. A method according to claim 12, characterised by 
10 determining distance to a single ground fault using a first 
order formula (27a, b,c} of the form: 

nnag[V_ AA _ p [3(7 M o ~ £qZabq)]* ) 
~ imag{(Z 1M / /W _p)[3(/ A10 - P 0 IyU»o)]*} 



16. A method according to claim 12, characterised by 
15 determining distance to a phase-to-phase ground fault using 

pre-fault measurements and a first order formula <28a,b,c) of 
the form: 

imag{VAA_ P Ml.AA0 -PqIabo)T ) 

0 )MLaao-£oLabo)T} 



20 17. A method according to claim 12, characterised by 

determining distance to a phase-to-phase ground fault avoiding 
pre-fault measurements and using a first order formula 
(29a, b,c) of the form: 

d= imaglQL, +V b )U.AAo- EqLabo)*1 

imag[Z llA (l a +Lb + 2/c 0 / /MO + 2k 0m LABo)(lMO ~£oI/iso)*l 
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18. A device for locating a fault from one end of a section of 
a power line (A-B) having means for receiving and storing 
measurements of current, voltage and angles between the phases 
at one first end (A) , means for receiving and storing a 
detection of a fault condition between said first and second 
ends <A,B), characterised by: 

-means for calculating symmetrical components of currents for 
said current and voltage measured at said first end, 
-means for calculating a distance (d) from said first end (2) 
to the fault (F) . 

19. A device according to claim 18, characterised by 

comprising: 

-means for determining a value for source impedance at said 
first end, 

-means for determining a value for source impedance at said 
second end. 

20. A device according to any of claims 18-19, characterised 
by comprising: 

-means for receiving a measurement of source impedance at said 
first end A. 

21. A device according to any of claims 18-20, characterised 
by comprising: 

-means for receiving a measurement of source impedance made at 
said second end B. 

22. A device according to any of claims 18-21, characterised 
by comprising means to receive a measured value (9) for remote 
source impedance at said second end (B) communicated by means 
of a communication channel (60). 
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23. Use of a fault locator device according to any of claims 
18-22, by a human operator to supervise a function in an 
electrical power system. 

5 24. Use of a fault locator device according to any of claims 
18-22, by means of a process running on one or more computers 
to supervise and/or control a function in an electrical power 
system. 

10 25. Use of a fault locator device according to any of claims 

18-22, to locate a distance to a fault in a power transmission 
or distribution system. 

26. Use of a device according to any of claims 18-22, for 
15 locating a fault on parallel power lines. 

27 . A computer program comprising computer code means and/or 
software code portions for making a computer or processor 
perform any of the steps of claims 1-17. 

20 

28. A computer program product according to claim 27 
. comprised in .one .or . more computer, readable *media. 

29. A data communication signal for locating a fault in a 
25 section of a power line included in a data transmission 

comprising a value of a measurement of source impedance made 
in respect of a remote and second <B) end of said section of a 
power line . 
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30. A graphic user interface for displaying a location of a 
fault in a section of a power line wherein a value is 
displayed for a distance (d) of said fault from a first end 
(A) of said power line. 



31. A graphic user interface according to claim 30, 
characterised in that the value displayed for the distance (d) 
is combined with a graphical representation of the relevant 
power line section or network. 

32. A graphic user interface according to claim 30, 
characterised in that the value displayed for the distance (d) 
is arranged to be displayed upon activation of a part of the 
graphical representation of the relevant power line section or 
network using a computer mouse or similar computer display 
selection means. 
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ABSTRACT 

A method to locate a fault from one end of a section of a 
power line (A-B) . Measurements of current, voltage and angles 
between the phases are made at a first A end of a said power 
line section. Upon detection of a fault condition between said 
first end and a second end of said power line, distance to the 
fault is found by calculating symmetrical components of 
currents for said current and voltage measure at said first 
end, then calculating a distance d from said first end 2 to 
the fault F the distance d to the fault using a quadratic 
equation of the form: 
B 2 d 2 + Bid+B Q -0 
where : 

B 2 = ^2_Rc^0O_lm ~ ^2_ImAjO_Re 

A value for source impedance at the first and/ or second end 
used in the method may be a representative value or a measured 
value. In other aspects of the invention a fault locator 
device for carrying out the method and a computer program for 
carrying out the method are described 
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